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1 About this note 

This is a note associated with my paper "From Quantum Dynamics to the Canonical 
Distribution - A Rigorous Derivation in Special Models" ( |cond-mat / 97072SBD . Here I 



describe all the technical details which are not discussed in the main paper. 

Please note that this is not (yet) written gular paper. I did not include any 

introductory materials or physical discussions. The proofs may not be optimally organized 
yet. 

The preset note is organized as follows. In Section ^ we prove a simple claim about 
the robustness of the non-resonance condition that we mentioned in the main paper. 
In Section ^ we prove the Theorem in the main paper. The theorem is essentially an 
application of the Chebysehv's inequality, and the proof is easy. In Section ^, we prove 
the Lemma in the main paper. As is clear from the table of contents, this is the hardest 
and the most technical part in our analysis. We have summarized the basic strategy in the 
beginning of the section. Section |^ is independent from the rigorous example we discuss 
in the main paper and in the present note. Here we deal with much more general class of 
models, and show that the "half" of the "hypothesis of equal weights for eigenstates" can 
be proved rather easily. 



2 Non-resonance condition 

Let us prove the statement about robustness of the non-resonance condition mentioned in 
the footnote [8] of the main paper. 

Let Ei be an eigenvalue of H and let the corresponding normalized eigenstate be 

n N 

^i^. = E E ^S;^)^. ® r.. (2.1) 
j=i k=i 

From the first order perturbation theory, we get 

dE,, 



dBy 



Ei^s;.')!'- (2-2) 



Suppose that the energy spectrum {Ei} violates the non-resonance condition for some 
EiS in the range Ei> en + 2A. More precisely, we assume that there are ii, Z2, u such 
that R = Ei^- Ei^ - {Ei^ - J = holds. 

We now shift all the Bk'^s in the lowest band (0, S) by a small amount, say d, keeping 
their spacing unchanged. If 

dR^dE^_m_dE^^dE, ^2 3) 

dd dd dd dd dd 
is nonvanishing, then we can conclude that the resonance is lifted for any small shift d. 
Also note that no new resonances are generated if we keep d sufficiently small. 

Let us assume that ^ happens to be vanishing. In such a (extremely rare) situation, 
we shift all the -B^/'s in the second band {6, 26) by d'. We can then say that 

dR_dEi_dE2_dE^ dEi 
dd/ ~ ~d(f ~ ~M ~ ~d(f ^ ~M 
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is nonvanishing. To see this we recall the representation ( |2.2| ). That If = means that 



there is a very special relation between J2 k' -Rawest J2]=i lv^(j^fc')P with = 1,2,3,4. Since 
(P{j,k) is determined as the solution of the Schrodinger equation (4) in the main paper, it 
has different decay properties in the "classically inaccessible regions" for different values 
of E. This means that Z]fe';second Sj=i \'^^{fk')\^ with /i = 1,2,3,4 cannot satisfy the same 
special relation as the corresponding quantities of the lowest band. So we conclude that 
the resonance is lifted by a small d'. 

The same argument works for the cases with multiple resonances. We see that all the 
resonances go away if we allow small d and d'. 

3 Proof of Theorem 

Let us prove the Theorem in the main paper. We first state and prove a general lemma, 
on which the desired theorem relies. Let be the eigenstate of the total Hamiltonian H 
with the eigenvalue E'. We assume that the initial state $(0) of the system is expanded 
as 

Ho) = Y.^E'^E', (3.1) 

E' 

and set 

7 = rnax Ije'I- (3.2) 

As in the main paper {A)^ denote the expectation value of the operator A of the subsystem 
in the state at time t. We have shown in the main paper that 



{{A),- (A),} <n\\\A\\J'f. (3.3) 

Then we have 

Lemma 1 Let A be an arbitrary operator of the subsystem. Let k > and A > be 
arbitrary constants. Then there exists a (k- dependent) constant T > 0, and we have 

where ta{T) is the total length of the intervals within <t <T at which 



(A), - {A),\ > A (3.5) 



holds. 



Proof: This is nothing but the Chebyshev's inequality, but we give a proof for completeness. 
For a function f(t) of t, we let 

1 



Mfit)] = -J^ dt{f{t)-f{t)f. (3.6) 
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Since ( |3l3| ) implies 



limt;^P)J<n^(Pl|jV 

1 Joo 



we see from continuity that for a given /s: > 0, there is T > such that 

VTU),]<{l + K)n\\\AWf. 

Now observe that 



X 



> A 



< 



(3.7) 



(3.^ 



(3.9) 



where the characteristic function x is defined by x[true] = 1 and x [false] = 0. By averaging 
roD over t such that < t < T, we find 



(3.10) 



which with ( ^.8| ) imphes the desired 



In order to prove the Theorem in the main paper, we have to evaluate 7 and choose 
appropriate A. 

We recall that we have the initial state of the form 



$(0) = ^^^^afclfc, 

k 



with ftfc nonvanishing only for k such that 

|i?-(en + 5.)|<y, 
for a fixed E such that + 3A < -E < i?max — 3A. To evaluate 7, we note that 

-iE' = ($i?',$(o)) = Y.w^(^k, 

k 



where we wrote 



n N 

j=i k=i 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



We shall prove in Section [4.6| that |v5(n,/c)| is less than 0(exp[— const. L^/^]) outside the 
interval {k^m, • • • , ^max} determined by the condition 



\E' - [en + Bj)\<X + const. L'^/l 



(3.15) 



For E' such that \E' — E\ > A + {en/2), the two ranges ( |3.12| ), (|3.15|) have no overlaps, 
and we see that Ije'I is small. 

So we assume \E' — E\ < A + {en/2). From ( |3.13| ), we have 

|7e'I < (max \ak\) V \ip{n,k)\- (3.16) 
k 



By using the above observation about |v5(n,fc)|, we find 

^max 

X!lv^{n,fc)l ^ H Iv5(n,fc)| + const. exp[-const.L^/^] 

1/2 



< I H \V{n,k)\ 



E 1 



const. exp[— const. L^''^] 



max ^min)^^^ + COUSt. exp [— COnst.L^^^] 



< j2\p{E' - e„ + A) + const, exp [-const. L^/^] 



By using the assume bound for at (see the main paper), we find that 

2c'A 



\lE' 



< 



(3.17) 



(3.18) 



and can set 7^ = 2c'A/e„. 

We recall that the desired theorem in the main paper does not specify the operator A. 
By linearity we can replace the phrase "for any operator A" in the Theorem by "for any 
of the operators A^^i, defined by {A^^i,)jj' = S^j6uj'" 

We now let A be one of the A^^^^^s, and apply Lemma |l] by setting 

/ \ \ 1/3 

A = n^PIL - , (3.19) 



and K = 1/2. Then we have 

T 



< 



A 



-2/3 



2ri2 \6r 



3c' fx 



71^ \ Er 



1/3 



(3.20) 



Since each of A^^j^ can have its own "bad" interval, the total length of the "bad" intervals 
of all A^^y is bounded by times the right-hand side of ( ^3.201 ). This is what appears in 
the right-hand side of (15) in the main paper. 

It only remains to estimate the systematic difference between the desired (^4)^^^° and 

{A)^. We first note that 



, can 
//3(F) 



^' rf/?(F) d 

E 



< 1 



dF {{AsY^^p^ - {AY^If) {^)7{f)) 



< 2^\E' ~ E\\\A\\^e, 



< 2||A||^7(^A + y )£„ 

< 2\\A\\^^{Sn)\ 



(3.21) 
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where 7 = dP{E + \)/dE. Thus we see 

<2P|L7(^nf. (3.22) 

By adding this new systematic error to the error a \\A\\^ in the (7) of the main paper, we 
finally get the statement of the Theorem. 



I A\ can I 1 2 / /t \ can 



4 Proof of Lemma 

In the present rather lengthy section, we shall prove the Lemma in the main paper, which is 
our main estimate. Throughout the present section, ci, C2, . . . , C23 denote positive constants 
which do not depend on L but may depend on {^j}. A, 5, and {Nr}. 

Let us give an outline of the present section. In Section [4.1| , we prove the statement of 
the Lemma, but assuming some new lemmas which will be proved in the latter sections. 



In Section |4.2| , we introduce the notion of "regular interval" , and study how the two 
index systems £ and (j, k) are related with each other. This is essential in getting the 
desired Boltzmann weights. 



In Section [4.3| , we decompose the whole region {1,2,..., nN} for the index £ into many 
subintervals. 

The next three sections are devoted to the estimate of the solution of the Schrodinger 



equation ( [4. 1]) in each of the above intervals. Therefore the topics of these three sections 
are asymptotic analyses of a discrete Schrodinger equation, and are not quite specific to 
the problem of deriving the canonical distribution from quantum dynamics. Although the 
techniques I use are not quite original, I present all the estimates since I could not find 



necessary estimates in the literature. In Section |4]J, we treat the solution in the "classically 
accessible region". The approximate solution is obtained by a speculation based on the 
quasi-classical analysis, and the difference from the true solution is rigorously controlled by 



the standard machinery of transfer matrices. In Section [4.5|, we treat the solution near the 



"turning points" where the quasi-classical analysis no longer works. We use the solution 
of (rescaled) continuous Schrodinger equation as an approximate solution, and control the 



difference from the true solution inductively. In Section O, we control the solution in the 
"classically inaccessible regions". 

In the next two sections, we extract information about the Boltzmann factor from the 



controlled approximate solutions of the Schrodinger equation. In Section |4.7| , we treat the 
"classically accessible region". We encounter an annoying phenomenon of "resonance", 
which locally inhibits the wave function to generate the desired Boltzmann factor. We 



suspect that this phenomenon is of essential character. In Section [4.8|, we get the Boltz- 



mann factor in the region where the wave length of the wave function is long. There are 
no resonances, and the proof is easy. 



In Section |4.9|, we fix some exponents, and complete the lengthy proof. 
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4.1 Proof 

We write the Schrodinger equation (4) in the main paper as 

V?£-i + ifi+i + 2ae,ipi> = 0, (4.1) 

with 

c^i = (4.2) 

where E is an eigenvalue such that 

en + 2X<E< S^ax - 2A. (4.3) 



We shall decompose the whole range of the index i into a disjoint union of Q intervals 

as 

{l,2,...,nN}= U 4, (4.4) 

1^=1 

where the intervals will be specified later in Section |4.3| . We here note that Ji and 
are special intervals which consist of the "classically inaccessible regions" with < —1 
and a£ > 1, respectively, plus small ranges in the "classically accessible region" attached 
to them. All the other intervals /2, /s, . . . , In-i are in the "classically accessible region" 
-I <ai < 1. 

We recall that we have two different index systems, namely, ^ with ^ = 1,2,..., nN, 
and {j,k) with j = l,2,...,n and k = 1,2, . . . , N, and these two are in one-to-one 
correspondence i ^ {j, k). We make the correspondence manifest by writing 

i^(j(i),kii)), or i(j,k)^ij,k). (4.5) 

This is a slight abuse of notation, but we hope there will be no confusions. 

The following is our main estimate. It sates that ipi produces the desired Boltzmann 
factor in each interval. 

Lemma 2 Let (pi he a normalized solution of ( j^. For each uo = 1,2, ... ,Q and each 
j = 1,2, . . . ,n, we have 



Ete/.x[j(^) =j] 



E 



< ciW^'^^L-^ (4.6) 



where the indicator function is defined by xitrue] = 1 and x[false] = 0. The exponent 
7] > will be determined later in Section (to be 1/12). Here 



is essentially the Boltzmann factor. The index i{uj) are taken from the interval I^^. For the 
intervals Ii and 1^, we take the corresponding £{1) and from the "classically accessible 
parts" of the intervals. 
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This lemma will be proved in subsequent sections by constructing local approximations 
for ip^. 

Given this lemma,we immediately get 
Lemma 3 Let ip^ he a normalized solution of For each j = 1,2, ... ,n, we have 



N 



k=l 



< CiWjL-\ 



with Wj satisfying 



and 



P{E - Bj - A) ^ ^ p{E - e, + A) 



E".=iP(i?-^,' + A) 



<w,< — 



piE - Sy - A) 



(4. 



(4.9) 



(4.10) 



Proof: This is almost trivial. By writing = J2eeiu, Iv^^P? have 

nN 



N 
k=l 



£=1 



UJ=1 



Since we have Y.2=i Pu; = ^ and P^ > 0, 



immediately implies ( [4. 81) with 



(4.11) 



UJ = 1 



(4.12) 



The normalization property ([4.9|) is trivial if we note ( [4. 71) and ( [4.12| ). The bound ( [4.10|) 
follows from ( ^ ) and (|4rT2| ) if we note E - X < t/^-^^^ < E + A for = 1, 2, . . . , and p 
is nondecreasing. ■ 



For any operator A on the Hilbert space T^s of the subsystem (with the matrix elements 
{A)jj' = (^j,y4^jv)), we define 



(4.13) 



We show that this expectation value is almost equal to the desired canonical expectation 
value 

,^.can _ Trs[Ae-/^^3] _ E-L,(A),,e-/^-^ 

^^f" ~ Trs[e-/3^^s] - Y.Ue-f^^ ' ^ ' 

This is an elementary estimate, and can be proved rather easily. 
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Lemma 4 For any operator A on Tis, we have 

(^)^-(^)r|ll^lloo{3/?7 + 7(^n)^} 
with 13 = (3{E) = d\ogp{E)/dE and 7 = d(3{E + X)/dE. 



(4.15) 



Proof: Let Ej = E — Ej. We shall prove upper and lower bounds for p{Ej ± X)/p{Ei) 
Since P{E) = dlogp{E)/dE, we have 



p{Ei) 



exp 



El 



dE(3{E) 



(4.16) 



for any E. By expanding 13{E) around Ei and recalling that d[3{E)/dE is increasing, we 
have 



< i (3{E) -{E- E,)(3{E^) < 1{E - E,] 

JEi Z 



(4.17) 



where 7 = P'{Ei + A). Substituting these bounds into ( |4.16|) , we have 

p{Ej + A) 



< e-/3(ej -ei)+/3A+(7/2) (sj -ei - A)^ 



< e 



-/3(ej-ei)+/3A+(7/2)(e„-£i)2 



and 



pjEj- A) ^ 
P(^i) " 



(4.18) 
(4.19) 



with p = p{Ei) for any j = l,2,...,n. By substituting (|4.18| ), ([4.19|) into (^4.10|) , we 
finally get 



-2/3A-(7/2)(£„-£i)2 



< 



< e 



2/3A+{7/2)(e„-£i)2 



(4.20) 



where Wj = e / Yl'j>=i e is the Boltzmann factor. 
We can finish the proof of (^4.15|) by observing that 



n 



J2iA),,iW,-W^ 



n 







> max 




1 ' 



< 



\A\\^{3pX + j{en-eif} 



(4.21) 



where the final line follows form ( [4.20| ) if 2/?A + (7/2)(£:„ — ei) < 0.7, which we shall 
assume. ■ 



To complete the proof of the Lemma in the main paper, we need one more estimate 



which will be proved in Section 4.6 
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Lemma 5 For any normalized ipi satisfying ( |^. we have 

N 

E WJ^^U',k) < const-e-^"-^^*-^''', (4.22) 
fe=i 

for any j ^ f . 

Given all the above estimates, to prove the desired Lemma in the main paper is straight- 
forward. With the expansion 

n N 

*s = EE<^0»*j ^Tfc, (4.23) 
j=i k=i 

the desired quantity becomes 

n f N ^ 

{^E, {A 1b)$e) = E (^)iJ' E W^VU',k) . (4.24) 
j,j'=i U=i ) 

Since the diagonal weigh Y.k=i lv^(i,fc)P is controlled by ([4.8|), and the off diagonal weight 
J2k=i ^ij,k)^{j',k) for J 7^ / by (H:.22|) , we immediately see that 

\{<i>E, (A® lB)$ir) - {A)^\ < 2ci \\A\\^L-^. (4.25) 

By combining this with the systematic error estimate ([4.15|) , we get the desired (7) of the 
main paper. 



4.2 Regular intervals 

We first recall special regular structure of the spectrum {Bk}k=i,...,N of Hb- In each interval 
((r — 1)5, r6) with r = 1, 2, . . . , i?, the energy eigenvalues Bk are spaced with exactly equal 
spacing br = 5{My.L)~^. This means that the whole index set is naturally decomposed into 
R intervals as {1, 2, ... , N} = Uf=i so that for any k E Kr, we have G ((r — 1)5, r5). 

Since Uq^^) = + B^, the structure of {Ui} inherits the above regularity of {-B^}- 
We say that an interval J C {1,2,..., nN} is regular if for any {j, k) ^ i E J, we have 
k G Kr(j) with ( J-dependent) r(l), r(2), . . . , r(n)(= 1, 2, ... , R). Thus, in the interval J, 
Ui is constructed by superposing n shifted copies of {Bk} each of them having exactly 
equal level spacings. 

The whole range of £ can be decomposed into a disjoint union as 

nR. 

{l,2,...,niV}= U J„ (4.26) 

s=l 

where each Jg is a maximal regular interval. We note that each regular interval has length 
of 0{L). We also remark that ( [4.26[ ) is not yet the decomposition 



We want to determine the behavior of the index £(j, k), Ui, and = (f/^ — E)/X in a 
fixed regular interval J, which is one of Ji, . . . , JnR- For simplicity, we write 

bj=brij), Mj=Mri,), (4.27) 
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for j = 1,2, . . . ,n. 

Let (j, k) i & J. Because of the regularity, we can write 



U(jM) = {k - Kj)bj + u 



(4.28) 



for some Kj and Uj (which are again J-dependent). The index ^ is determined by ordering 
U{j^k) so that Ui < Ui+i. 

To get an exphcit formula for k), we count the number of (j', k') G J such that 



For a fixed j' ^ j, the number of k' with ( |4.29| ) is 



(4.29) 



(4.30) 



where [■ ■ ■] is the Gauss symbol. Summing up these contribution as well as that from the 
indexes (j, k') G J, we get 



i{j,k) = i,-l + {k-K,) + l + Y^ 



i'=i L 



{k - Kj)hj + Uj - Uj: 



(4.31) 



where we used ( [4.28|) _, and is the smallest index in J. For the latter use, we substitute 
the relation hj = 6{MjL)-^ into (|OTl ) to get 



k + r]j,. 



(4.32) 



where £q and rjj j/ are constants (which may depend on L). 
From ( [4. 311 ), we see that 



with 



Let 



i{j,k)-{n-l)<i{j,k)<i{j,k), 



j'=i 



Note that 



L 



-1 9^ 



(4.33) 
(4.34) 

(4.35) 

(4.36) 
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where 



9 = ^ (4.37) 



is an L-independent quantity. Observe that 



Uj — u 



= U^j,k)-u, (4.38) 



where 

u = bY,^ (4.39) 

is an energy near the bottom of the interval J. From ( [4.38| ) and ( [4.33| ), we get 

U + U <Ue<b£ + U + b{n-l), (4.40) 

for i in the regular interval J, where U = u — bio. 
We introduce a linearized by 

U + U-E i U-E 

= ^ =^Z + ^' (4.41) 

for I E J . Then from (|4.40|) , we find 

I _ I . b{n-l) g{n-l) 

\oii-a{,\< = . (4.42) 



4.3 Decomposition into intervals 



We describe precise definitions of the decomposition ( |4.4| ) of the intervals. The intervals 
Ji, . . . , /q are properly ordered, and covers the whole range {1, 2, . . . , nN} without any 
overlaps. 

We define the first turning point £t as the minimum £ such that > —1. Then the 
first interval is defined as 

/i = {l,2,...,£t+[c2i^'/']}. (4.43) 
Again [■ ■ ■] is the Gauss symbol. Similarly we define 

In = {i[-[c2L^'%....nN], (4.44) 

where the second turning point is the maximum I such that < 1. 

The intervals I2, ■ ■ ■ , Ir all have the length = [csL^^^^^^"^]. T is determined as the 
minimum number such that Y,^=2 ^ [04^^"^]. The exponent e > will be determined 
later, but it must satisfy 

l-e>^-r], (4.45) 

because we must have I]5'=2 \^^'\ — \^^\- Similarly the intervals /^-r+i, • • • , -^f^-i have the 
length [c3L(^/3)-»?]_ 
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The remaining intervals /r+i, • • • , Iq-f are defined to cover the (wide) remaining region 
[it + [02^1/3] + . . . , A - [C2L^^^] - [caL^'']}. We require for = T + 1, . . . , ^] - T 

that has the length 

[c,L'-^'/2] < m < [c,L'-'% (4.46) 

where 6 is an exponent to be determined later, and is contained in a single regular 
interval defined in Section [4.2| . These two conditions are easily satisfied since the lengths 
of the regular intervals are of 0{L). 



4.4 Approximate solutions in the "classically accessible region" 



We will construct approximate solutions of the Schrodinger equation ( ^4.1|) in the intervals 
I2, ■ ■ ■ , In-i, which are within the "classically accessible region". 

For simplicity, we denote by / = {£1, £1 + 1, ■ ■ ■ , £2} one of the intervals I2, ■ ■ ■ , In-i- 
Since the interval I is entirely contained in a single regular interval J, we have correspond- 
ing ae defined by ( [4.41| ), which is a linear approximation to ae. We also write 

a = ai>2, f3 = ae^. (4.47) 

We start from an abstract theory for (rigorously) evaluating the difference between an 
approximate solution and the true solution of (|4.1|). Let ipe be an approximate solution of 



(4.1) with a£ replaced by its linearlization ae in the sense that 

ipe+i + ipe.i + 2aiipi = Se (4.48) 

holds for i & I = {£1 + 1, ... ,£2 — Here Si is an error term, which must be small. 

For a solution ifg of the original Schrodinger equation ( |4.1| ), we write the deviation 
from the approximate solution as 

fi = ipi- ^pe- (4.49) 

From (|0| ) and ( [4.48| ), we find 

+ + 2aefe = a,, (4.50) 

with 

(Ti = 2{ai - ae)i!e - 5e. (4.51) 
As usual the equation (|4.50| ) can be written in a matrix form as 

where the transfer matrix is 

J,= {-f' -;). (4.53) 

Considering the second order nature of the equations ( [4.1| ) and ( [4.4(^ ), we can assume 
without loosing generality that /^^ = fe^^+i = 0. Then by inserting ( |4.52| ), we find 

E T,T,_,...T,^,f"^^). (4.54) 



k=ei+i 
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Since we have \ai\ < 1 in the "classically accessible region", the transfer matrix Ti of 
( [4.53| ) has two eigenvalues e£,e£ with \e(\ = \ei\ = 1, where 



ee = -Oi + i\Jl- ae"^. 
Thus there exists a regular matrix such that 





Let us define 



(M,)-iT,M, 







1 
1 



(4.55) 



(4.56) 



(4.57) 



which is expected to be small since and M£_i are very similar with each other. Then 
from ( |4.56| ), we have 

T^T^-i ■ ■ ■ Tfc+i 

^M,(;' °)(M.,-'M„(-- ^)(M,_0-'- 

■••(l + Dt«)(Y JL)(M»+,)-'. (4.58) 

For any a,b = 1, 2, we denote by (A)^^;, the a, 6-component of a 2 x 2 matrix A. We 
now assume that 

\{Di)aA < di, (4.59) 
for any a,b = 1,2. Then ( [4.5^ ) implies 



n (1 + 2d, 



j=k+i 

exp 



|(T,T,_i---Tfc+i),,fe| < ||M,||^ (Mfe+i)-i 

< ||m,il||(m,+o-' 

To make the estimate ( 4.60| ) more concrete, we fix precise form of as 



2 Y: d, 

j=k+i 



(4.60) 



whose inverse is 



(M,)-i 



2^2(1 -a/)i/4 ^-1 ei 



(4.61) 



(4.62) 
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Since defined by ( [4.57| ) is vanishing if ai = a^-i, we can evaluate by expanding 
it in («£ — a;£_i). Tlien we find tliat essential contribution comes from the first order in 
the expansion, and the desired bound ( [4.59| ) is satisfied with 



df 



provided that 



I- afP' 
2 



1 - a/ > 



(4.63) 



(4.64) 



which is automatically satisfied in the present interval. Substituting ( [4.63|) into ( |4.60|) , we 
find 

2{ai - ak+i) 



O'i^e-i ■ ■ ■ Tfc+i)a,fe| < II M| 



'M 



k+l) 



v-1 



exp 



1-72 

for any a,b = 1, 2, where 7 is one of with i E I which gives the smallest 1 — ag^. 
At this stage, we impose a condition on / that 

a — P 



(4.65) 



1 -Y 



< C6, 



(4.66) 



where a,P are defined by ( |4.47] ). Then ( |4.65| ) simplifies as 



(T^T^.i ■ ■ ■ Jk+l)a,b\ < 



IM, 



M 



=C6 



< cr(l - «/)-i/^(l - a^)-'/\ (4.67) 
where we used (|]6l|) and (g^H). Substituting (^M ) to (|43^ ), we get a useful bound 

e 

C7 



e+i 



< 



(l-a,2)i/4(i_y)i/4^E^J^'=l 



(4.68) 



This completes the general theory. 

We now explicitly construct an approximate solution ijj£ which appears in ( |4.48| ). We 
first define d by 

cos Q = —ag. (4.69) 

Since ag is increasing in i and satisfies |a^| < 1, we see that Q is increasing and satisfies 
< < TT. Let us define 



£-1 

^£'=£1 



2 



Then we can write down our approximate solution as 

cos(Zf + 



1p£ 



(4.70) 



(4.71) 



where ^ is an arbitrary constant. The form ( [4.71|) is (heuristically) obtained by following 
the standard idea of quasi-classical analysis. We learned that this approximate solution 
was written down long time ago by Bethe^. 

2 H. Bethe, Phys. Rev. 54, 955 (1938). 
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We now need an estimate for the error term 6i which appears in ( [4.48| ). To do this, 
we substitute the concrete form ( [4.71|) into the left-hand side of ([4.48|) , and expand the 
resulting quantity in a power series of (a^+i — ae) = (a^ — ai^i) = b/\(= g/L). (See 
( [4.41| .) We find that the first order terms cancel out, and the essential contribution comes 
from the second order. We skip the tedious but straightforward calculation, and describe 
only the final result, which is 



< 



Oil 



2(1 - (a,)2)3/4 (1 _ («,)2)3/2 



(a,)2)5/4 ^A(l-(a,)2)i/2^ 

/T\ 2 



< 



(1-(«,)T/' \\, 



(4.72) 



Here we used the fact that ~ in the sense of ( [4.42|) . 

By recalling the definition ( |4.51| ) of cr^, and using ( [4.72| ) and ( |4.42| ), we find 



< 



2n 



;i-«,2)i/4 ^a; 

where we noted that sin^^ = \/l — a^^. Since 



+ 



^ b 
k=ei+i ^ 



we find 



El I ^ V 



2n{a-p) 3{a-/3) (b 



+ 



Going back to (|4.68|), we finally see that the relative error is bounded as 



(4.73) 



(4.74) 



(4.75) 



< 



C7 



/ ) fc=£i+l 



< 



2nc-j{a — [3) ?>c-j{a — (3) ( b^ 

(1-72)1/2 + (1 _ ^2)5/2 [^^ 

For the latter uses, we want to get a bound of the form 



fiJ sin Ci 



< CsL-\ 



(4.76) 



(4.77) 



with a constant cs which does not depend on L and the specific interval. 

We shall write down conditions that the exponents t], e and 6 should satisfy to get 
( [4.77| ). We start from the case where / is one of the second type intervals, i.e., /^^ with 
uj = 2,...,roTUJ = Q — r + l,...,Q — l. These are the intervals which are relatively close 
to the turning points. For such an interval, we have 1 — 7^ > C2L^^^{b/ X) = C2gL~'^^^. 
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(See ( |4.37| ) for the definition of g, which is an L-independent quantity.) We also find 
[a - 13) < C3L(i/3)-''(6/A) = C3^L-(2/3)-r?_ It is easy to check that the condition j WM) is 
satisfied. By substituting these bound into ( |4.7(j| ), we see that 



/^VsinO 



(4.78) 



Thus we find that the desired bound ( [4.77| ) is indeed satisfied. 

We then consider the case where I is one of the third type intervals, i.e., with 
= r + 1, . . . , - r. Here we have (1 - 7^) > CiL^~^{h/\) = c^gL'^, and (a - p) < 
c^L^~'^^ (b / X) = c^gL""^^. Therefore the condition ( [4.66| ) is satisfied if 



-2e + £<0. 
Substituting these bounds into ([4.76|) , we get 



(4.79) 



fey sin Ci 



V C4 vc4^g 
3C7C5 



C4 Vc4^g 



_^2e-l 2,-2e+(e/2) 



Therefore, we have the desired bound ([4.76 ) provided that 



and 



1 



29 + - < -71. 
2 - ' 



(4.80) 



(4.81) 



(4.82) 



Later in Section |4.9| , we determine all the exponents appear in the proof. As these expo- 
nents, we will set ?7 = = 1/6 and e = 1/3, which satisfy ( |4.81D and ( [4.82[ ). 
To summarize, we have proved 

Lemma 6 Let ipi he a real solution of ( |^. j[ ). Let I he one of the intervals with uj = 
2, . . . ,Q — 1. There are real constants A and ^, and we have 



-A 



cos{Zi + 



<C8 



A 

A/sin (i 



(4.83) 



for £ e I. 



4.5 Approximate solutions near the "turning points" 



We will construct approximate solutions of the Schrodinger equation ([4.1|) in the "classi- 
cally accessible parts" of the intervals Ji and Jq. The remaining "classically inaccessible 
parts" will be discussed in Section |4.6| . Near the turning points, the wave length of the os- 
cillation of ife becomes long and the quasi-classical approximation becomes useless. Instead 
we try to approximate ipe by the solution of a rescaled continuous Schrodinger equation 



that corresponds to (]4J 
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Let us discuss the interval Ji which contains the first turning point £t at which ~ 1. 
The treatment of Jq is essentially the same. We want to construct approximate solution 
of (fii for i in the interval Ji = j^t, 4 + 1, • • • , 4 + [c2-^^^^]} C h. let us write 



= 2{ai + 1 

and rewrite the Schrodinger equation ([4.1| ) as 

ipe-i - 2(pi + (pi+i + 
From (|CT| ) and j ^^) , we find that 



2(n-l)b , , 1 
< ^ ^ ^ =2{n-\)gL-\ 



(4.84) 



(4.85) 



(4.86) 



We want to approximate by a continuous function ?/'(a;) with x = L ^^'^(i ~ 
With this correspondence in mind, we divide (H:.85|) by (L~^/^)^ to (roughly) get 



We then note that 



L-2/3/5,^L-2/3|(£-4) = ^x. 



(4.87) 
(4.88) 



So we are motivated to consider the continuous Schrodinger equation 

^"{x) + gx^{x) = (4.89) 

Let %l){x) be the solution of ( [4.89 ) in the region x > 0. By an explicit calculation, one 
finds that two independent complex solutions of ( [4.89| ) are given by 



^{x) = ^h[{2/d,)^x^'\ 



and ip{x), where 



h{z) = H^^l{z) = J,„{z)+iY, 



1/3 l^J 



(4.90) 



(4.91) 



is the Hankel function (or Bessel's function of the third kind). From the asymptotic 
behavior of the Hankel function, we find 



i){x) ^ x~^/Sxp [i {(2/3)^/^x^/2 - (57r/12)} 



(4.92) 



for X ^ 1. 

From now on, we use the index m = i—it for convenience. To control the approximation 
rigorously, we first Taylor expand ip(L~^^^{m ± 1)) to get 

V^(L-i/3(m + 1)) - 2'ilj{L-^/^m) + i>{L-^/^{m - 1)) = L-^/^ip"{L-^/^m) + L'V^, (4.93) 

where 



6 



(4.94) 



18 



with L-^/^{m + 1) < x' < L-^/^m < x" < L-^/^{m - 1). By using and (|]93|), we 

find 

^(L-^/3(m + 1)) - 2^lj{L-^'^m) + ^lj{L~^''\m - 1)) = -\m^{L-^'^m) + L'V^, (4.95) 

A 

where we noted that L~'^/^ gL'^^^m = (b/ \)m. We divide the equation ( [4.95| ) by il'{L^^^^{m+ 
1)) to get 



1-2 



+ 



V^(L-V3(rn+ 1)) V^(L-i/3(rri+ 1)) 
a"^^(L-V3(^ + 1)) + 



V'(L-i/3(m+ !))• 



(4.96) 



We now let iIj{x) be the specific complex solution ( |4.90|) , and try to control a complex 
solution of (|4.1|) such that ^it+m — ^{L~^^^m). Information about the desired real solution 
can be read off easily, as we do at the end of the present section. 

Let us introduce a complex quantity Fm by 



Then the Schrodinger equation (|4.1| ) becomes 



(4.97) 



F^+ML-'^Hm + 1)) - 2F^tlj{L-'/^m) + F,r.-ML~"\m - 1)) 
= -l3i,+mF^^{L'^'''m). 

We divide this by FmipiL^^^^im + l))to get 

Fm+i ^ tPiL-'/^m) ^ ^{L'^^im - 1)) 



(4.98) 



-(3i 



^(L-i/3(rn+ 1)) F„ tlj{L-^/^{m + l)) 



XL-V3(rn+ 1))' 
From ( [4.96 ) and ( 4.99| ), we get the following recursion equation for F„ 



(4.99) 



F, 



m+l 



- 1 



1 



p 



A ; ^(L-i/3(rri+ 1)) ^/>(L-V3(^ + i)) " 
We now use the asymptotic behavior ( [4.92| ) of ipi^x) to see that 



(4.100) 



i){L-^l^{m-l)) 



^(L-i/3(m + 1)) 



< 1 + const. 



ip\x) 



^-1/3 



iIj{x) 

< 1 + (Cg + Cioxl/2)^-l/3 

= I + C9L-I/3 + C10L-I/W/2, 



(4.101) 
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m] 



V^(L-V3(m + 1)) 



<l + CoL-V3 + CioL-l/2^1/2 



(4.102) 



and 



^(L-i/3(^ + 1)) 



< const. 



Let us define 



< cii + ci2a;^/^ 

= cii + ci2L-i/W/^ 



m—1 



(4.103) 



(4.104) 



By using the estimates ( |4.101| ), (|4.102| ), (|4.103|) as well as ( [4.86|) , the recursion relation 
( |4.10q ) reduces to 



\G. 



■m+l I 



< 



\Gr. 



\G.n 



where we used 



m— 1 



- 1 



(G„ + l 



1-1 



(4.105) 



(4.106) 



We now assume that Gi = 0. Then we can estimate \Gm\ by repeatedly using ( [4.105[ ). 
Let us assume \Gm\ < G holds for any m with 1 < m < C2lv^/^ where the constant G will 
be determined later. Then ( ^.105| ) implies for any m with 1 < m < C2-L^/^ that 



iG^i < T.{m-\Gj-i\} 

< G(C2C9 + C2^/'cio) 

+G^(c2L^/3 + C2C9 + C2^/^Cio) 
+C2Ci3(L-2/3 ^ C2C9L-I + C2^/'cioL-l) 

+C2C11 

L-2/3 + c2^/2^i2L-2/3. (4.107) 
G = Cl4L-2/^ (4.108) 



We now set 



and substitute this relation into ( |4.107| ). We then find that, for sufficiently small (but 
L- independent) C2, ( [4.107| ) reproduces \Gm\ < G with the same G. This proves the upper 



bound 

for m with 1 < m < 02^^/^. 



\Gm\ < C14L 



-2/3 



(4.109) 
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Assuming Fq = 1, we finally get 

IF - II 



' F 



< expfmcuL"^/^] - 1 



(4.110) 



for m with 1 < m < 02-^^/^. Recalling ( [4.97 ), we have established the desired relation 



To control the desired real solution of (|4.1| ), we only have to sum up the complex 
solution and its complex conjugate with appropriate complex weights. This proves 

Lemma 7 Let ipi be a real solution of Then there is a complex constant A, and we 

have 



<fe - {A^{L-'/-\i - A)) + A^{L-^l\l - 4))} 



(4.111) 



for £ G Ji = l^t, • • • , A + [c2-^"'^^'^]}; where ■\\}{x) is explicitly given by ( \i-9(] ) 



We still have to treat the solution in the "classically accessible part" of the interval In, 
i.e., Iq = |£[ — [02^^'^^], . . . Since the analysis is exactly the same as that for Ji, we 
only present the final result. 

Lemma 8 Let (pi be a real solution of U-A )- Then there is a complex constant B, and we 
have 



- {-lY {i?^(L-i/3(4 - £)) + i?^(L-V3(£^ _ 
<2c,,\B\ ^{L-'/'{i[-e)) 



(4.112) 



for i E In, where ip{x) is explicitly given by ( U-dOj ). 

4.6 Decay of the solution in the "classically inaccessible regions" 

We now study the solution in the "classically inaccessible region", which is characterized 
by \ai\ > 1. Since the Schrodinger equation ( |4.1| ) imphes 



2a, 



we get a convexity inequality 





+ 




2 







< 



(4.113) 



(4.114) 
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This means that \(pi\ cannot take a local maximum. 

Let us focus on the region = + 1, . . . , nN} C Iq. Then ([4.114|) means that \ip(\ 
is decreasing because of the boundary condition (fnN+i = 0. Then ( ^.1141 ) further implies 



< 



\ai\ 



and hence 



We note that ( [4.42| ) imphes 



(4.115) 



(4.116) 



ng 



<exp [-cie{9/L){i'-i[-n)], 



(4.117) 



for i' — i[ < cijL^ ^ and any /i > 0. Recall that g is L-independent as in ( |4.37| ). Substi- 
tuting (|TTl3) into (|4.116|) , we get 



< l^j^Jcisexp 



0^ 



(4.118) 



for i' — i[< Ci-jL} ^. This means that \ip^\ decays very rapidly in the "classically inacces- 
sible region". For ^ > cijL}~^, we have 



< Ci8 exp 



(4.119) 



In the most applications (see (p.l5|) and Lemma 13), we set /x = 1/3 in ([4.119|). 



4.7 Boltzmann factor from the interior of the "classically acces- 
sible region" 

We are now ready to prove the most important Lemma 0. We first treat the intervals 
with cij = r + 1, . . . , n — r. These intervals are located in the interior of the "classically 
accessible region", where the wave length of ipi is relatively short. In such situations, we 
must face possible "resonances" between the oscillation of \ipe\'^ and the quasi periodic 
behavior of £(j, k) (with j fixed and k varied), which (locally) destroys the desired "equal 
weighted" behavior. We will prove that such resonances are located in short intervals and 
do not have significant contributions. 

Let I be one of 1^) with C(j = r + l,...,fi — F. Our final goal is to evaluate the quantity 
Sj/iJ2]'=i Sj') for j = 1, 2, . . . , n, where 

S, = j:x[j{i)=j]\v^- (4.120) 

Because of the definition of the interval, we have 

l-a/>C4^L-^ (4.121) 
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for any £ G /. (See the end of Section [4.4| .) We also recall that the length of I satisfies 

< |/| < [csL^-'^ (4.122) 
To evaluate the sum Sj explicitly, we further decompose I into subintervals as 



^ = U 

q=l 



(4.123) 



with each Iq having the length |Jg| = [cigL''], where z/ > is another exponent to be 
determined later (to be 1/3). By iq we denote the smallest element in Iq. Consequently, 
we write 



Q 

E 

q=l 



with 



= T.XW)=3]\^ 



(4.124) 
(4.125) 



Let i & Iq. We note that 



sinO 



sin 



1 - ai 



< max(l — a^^) ^ 

~ C4 

< ^L-^ 

C4 



where in the final line we used a new assumption on the exponents 

u + e — 1 < —7]. 



(4.126) 



(4.127) 



We now use the estimate ( [4.126| ) and the approximate solution ( [4.83|) to write the sum 
( [4.125| ) more explicitly as 



(a) 



A' 



sinQ 



(4.128) 



1 lain 



where R^^'' satisfies 



R 



(<?) 



— • /■ ' 



(4.129) 



where /^■'^ is the subset of Iq with j(£) = j, and |/^-'''| denotes the number of its elements. 
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We now want to evaluate the sum over cos^ terms in ( [4.128| ). Let Cg = Ci ) 



Zq = Then we have 



:cos[Z, + e])' - (cos[Z, + e + - iM 



< 2 max 



Ze-{z, + C,ii-i,)} 

< 2 max 10 - Cq\ \Iq\ 

\ ) 

< 2 (max(l - ae'^y^^A 



- as 



< 2 j^(ci9)'L2'^+(^/2)-i 

V C4 

< 2. &cl9)2L-^ 

V C4 

where the final line again makes use of the new assumption 



(4.130) 



2i^ + - - 1< r?. 
2 - ' 



(4.131) 



We introduce a new constant C,q = Zg + C, — (qiq- Then ( [4.1301 ) essentially means that 
(cos[Z^ + ^])^ — (cos[Cg^ + ^g])^- So we are motivated study the sum 

T.X[j{i)=j]{cOs[Cqi + Qf 



+ 1 ^ = j] (e^^iCq^+'^q) + e-2^(f^^+€«)) , (4.132) 



which is a good approximation to the sum in ( |4.128| ). We expect the sum over oscillating 
exponential in ( |4.132| ) to be small, but this is not straightforward. Indeed, the sum is not 
at all small if a "resonance" between x[j(^) = j] ^ind the exponential term takes place. By 
using ([4.32|) , we rewrite the sum of the first exponential term (times e"^*^'') as 



^ exp 2</; + 2tQ 
ks.t.e{j,k)eiq \ 

^max ^min / 

J2 exp 2</; + 2<g E 

p=0 \ j'=l 



k + 7]jj 



(4.133) 



where [■ ■ ■] is the Gauss symbol. The sum in the second line is over k such that £(j. A;) G Iq 
for the fixed j, and we denote this range of k as \kmm-, ■ ■ ■ , ^max}- We let s = [(fcmax — 
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^min)/Afj], and write p = Mjs + r. Then the above sum becomes 



s=0 r=0 \ j'=l 



+e- 



E 2zCg E 



p=sMj 



j'=l 



Mi 



P + Vk 



where 



. s-lM,-l / ri ( 

s=Q r=0 \ j'=l [ 

+ (the same second term) 

- - 1 _ p'iiC.qMs ( n 

+ (the same second term), 



Mi 



M^ 
Mi 



(4.134) 
(4.135) 



By taking the absolute value of ([4.133|) and (|4.134| ), we find 



Ex[j(^)=j> 



< 



< 



sin^^Ms| 
|sinCgM| 

2M 



1 Mi 



sinCgMl 



(4.136) 



As we have anticipated, the right-hand side of (|4.136|) is usually small (compared with 
|/^-'^|), but becomes large near the "resonance" points where (qM is equal to an integer 
multiple of 27r. 

To control the final sum ( ^.124| ), we classify the sub intervals Ii,...,Iq into "good" 
ones and "bad" ones. A subinterval Ig is said to be good if 



|sinC,M| 



(4.137) 



and to be bad otherwise. We want to know the possible number of the bad subintervals. 
First we note that the perfect resonance QM = 27r x (integer) can take place in the whole 
interval I at most once, since Q varies at most by 0{L'^~^^^^'^^~^) within /. We suppose that 
there happens to be a perfect resonance (with (^res) within J, and see how many subintervals 
around it are "infected" and become bad. From ( [4.137|) and that \Iq\ = [cigL'^], we get 



|C<J Crcsl ^ 



res I ^ ; 
C20C19 



(4.13^ 
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for C^q corresponding to a good Iq. This in turn means that the total number of the bad 
subintervals is bounded from above by 



C20C19 V<?/'G/ 



max|0-0,| Q < {gc,L~''r'Q 



C20C19 



< 



1 

fi'C5C20Ci9 



where in the final line we assumed that 

7] — i> + 26 < —1]. 



(4.139) 



(4.140) 



To summarize, we have proved the following for the desired sum (im. When 



Iq is a "good" subinterval, we have 



S. 



(<?) 



A' 



2 sin Q 



< C21- 



sin 6 



-\mL-^. 



(4.141) 



When /„ is a "bad" subinterval, we have essentially no control on 5',-'^\ and can only say 



(9) 



< C22- 



A' 



sinO 



|/(i)| 



(4.142) 



but we have the estimate ( [4.139| ) for the possible number of the "bad" intervals. Recalling 
(|]12|), we sum up (|4.141| ) and ( |4.142|) with (|4.139|) in mind. We then get 



< C23S,L-\ 



with 



This immediately leads us to our goal (|4.6|) that 

— ^ _ < c,W,L-\ 
2^j'=l ^j' 

with 



(4.143) 
(4.144) 



\tU)\ 

Kg I 



(4.145) 



(4.146) 



which is independent of q and is equal to ([4.7]). 
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4.8 Boltzmann factor from regions with long wave length 

We still have to prove the main Lemma ^ for the interval I^^ with uj = 1, . . . ,r and 
u = Q — r + 1, . . . ,Q. In these intervals, the wave length of (pe is quite large, and the 
resonance effect (which made the estimates in Section |4.7| difficult) does not take place. 
The proof is rather straightforward, and we will be sketchy here. 

In these intervals, we use one of the Lemmas |^, |^ or ^ to get controlled approximations 
for the solution ipi of the Schrodinger equation ([4.1| ). In most of the situations, \(pe\ can 
be treated essentially as a constant, and the estimate of ( [4.12(j| ) becomes as trivial as 

= T.xm=j]\Ve\' ^ \Vi\'T.x[m=j]. (4.147) 

eel eel 

The worst case that we have to worry about is when if£ changes its sign inside /. To 
investigate such a situation, we approximate ipe ~ {i — imin), and consider a small interval 
I = {^min, • • • , ^max}- The sum to bc evaluated is 

Sj= E X[j{^)=m-^minY. (4.148) 

Though may be a very complicated function of i, ( [4.32|) implies that it is (at least) 
a periodic function of i with the period M = jyj=i Mj. Assuming £max — ^min = sM, we 
evaluate for r < M the sum 

s=0 



to see that the (unwanted) r-dependent terms are of order 0(l/s). Thus in order for 
Sj/ J2f=i Sj' to become the desired Boltzmann factor with relative error 0{L~^), we must 
have that s > 0{L^), and hence 

Cax - Cin > const. L''. (4.150) 

In the present situation, ^^ax — ^min is determined by (the smaller of) the length of the 
interval and the wave length of the oscillation of ipi. Since the wave length in this region 
is longer than a constant times L^/^, the required conditions for the exponents are 

1>V, l-V>V- (4.151) 
These guarantee the Lemma for the desired intervals. 

4.9 Determination of the exponents 

In order to complete the lengthy proof of Lemma 0, we have to fix the values of the 
exponents rj, e, 9 and v so that several requirements are satisfied. 



g(2s + l)(s + l) 
6 

-s{2s + l){s + l)^~. 



+ s{s + l)Mr + {s + l)r^ 
M\l + 0{l/s)), 



(4.149) 
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We now recall the requirements about the exponents appeared as ( [4.45| ), ( [4.79| ), ( [4.45| ) 
( TO ), (PI, (im), (mil), ( ^140| ), and (liTsTI) , which are 



e 

-20 + - 
2 

p + e-1 
2 





1 


> 


3 


< 


0, 




1 


< 






2 


< 


— 71 

'h 


< 


-V, 


< 


-V, 


< 


-V, 


> 


V, 


> 


V- 



e 
2 

l-V > V- (4.152) 



As a solution which satisfies all of ([4.152|) , we shall choose 



This completes the proof. 



5 General upper bounds for \^(j^k)\ 

In the main paper, we noted that the "hypothesis of equal weights for eigenstates" can be 
partially proved. Let me describe precise statements. 

Recall that the "hypothesis of equal weights for eigenstates" consists of two parts; 

1. |v5{j,fc)P is negligible in the "classically inaccessible region." 

2. |v5(j,fc)P takes appreciable values all over in the "classically accessible region", and 
the value is essentially determined by a function / of the energy E — {ej + Bk). 

Clearly the second point is much more subtle. In fact if we take models with certain 
conservation laws, the second point is easily violated. The main point of our (unproven) 
hypothesis is that the above 2 holds in a general model which do not have any special 
symmetries or conservation laws. 

On the other hand the first point is more universal, and may be stated for a large 
class of models rather easily. Here we prove two results which justify this second point. 
There can be various theorems under different assumptions. We here present two simple 
theorems, which are in some sense complimentary with each other. 

We use the same notation as in the main paper, but now take general coupling Hamil- 
tonian H', and denote its matrix elements as 

Vj,k;j',k' = (^i ® Tfc, H'^^j, ® Tk') . (5.1) 
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We also denote by Hs and Hb the Hilbert spaces for the subsystem and the bath. 
As in the main paper, the the Schrodinger equation is 

{Hs ® 1b + Is ® i^B + H')^ = E^. (5.2) 

Again expanding the eigenstate as 

n N 

= EE<^o»*j®rfc, (5.3) 
j=i k=i 

^ is rewritten as 



{E - {sj + Ek)} (p(_j^k) = E ^j,k;j',k' V^{j',k'), (5.4) 

j',k' 

for any j = 1, . . . ,n and k = 1, . . . , N. 
We define 

A = sup ^ |V,-fc;j/,fc/| . (5.5) 

j,k 

The first result is very simple. 

Theorem 9 Let a normalized state $ G Hs^Ti-B satisfy ( \5. 2J . Then the coefficients (p(j,k) 
defined by ( \5. 3i ) satisfy 



Proof: From (|5.4| ), we see that 

\'PU,k)\ < 





k' 


Vj,k;j',k' 


\fij',k')\ 




E 


- {ej+Bk)\ 



where we noted |v5(j,fc)| < 1 because $ is normalized. 



The bound ( |5.7| ) is very crude, but gives us the idea that |v5(j,fc)| may be large for 



\E — (sj + Bk)\ ^ A, and is small for \E — {sj + Bk) \ 3> A. 

To get stronger result, we further assume that there is a constant D > such that 
Vj^k;j',k' = if \{ej + Bk) — {6jr + Bk')\ > D. In other words, the interaction Hamiltonian has 
no matrix elements between the basis states whose unperturbed energies differ by more 
than D. The assumption may hold for some interaction Hamiltonians which represent 
(near) elastic scattering, but not for some interactions. Under this assumption, we have 
the following. 



Theorem 10 Let a normalized state $ G Tis ® 'Hb satisfy (\5. Sj ). Then the coefficients 
'^{j,k) defined by ( |3T^ ^ satisfy 



^ \ ri{j,k) 
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where the integer n{j, k) is defined as 



n{j, k) 



\E 



D 



(5.9) 



Here [■ ■ ■] is the Gauss symbol. 



Proof: The bound ( |5.8|) (with the convention 0! = 1) is trivial for j, k such that n(j, k) = 0. 
Assume ( ^.8|) for all j', k' such that n{j',k') < n — 1. Take j, k such that n{j,k) = n. 
Then from (14. 11) we have 



< 



< 



A 



max 



\E - (Si + Bk) \ ij',k')s.t.\ie,+Bk)~ie,,+B^,)\<D 



\V{j',k')\ 



X 1 



A 



(n-l) 



nD {n - 1)! \D ^ 



< - - 

- n\\D. 



(5.10) 



which proves the desired bound. 



The bound ( |5.8| ) shows that actually decays very rapidly in the "classically 

inaccessible region." 
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